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Abstract. We examine a large class of inhomogeneous spherically symmetric spacetimes that gen- 
eralize the Lemaitre-Tolman-Bondi dust solutions to nonzero pressure ("LTB spacetimes"). Local 
covariant LTB objects can be expressed as perturbations of covariant quasi-local (QL) scalars that 
satisfy evolution equations of equivalent Friedman-Lemaitre-Robertson-Walker (FLRW) scalars. 
Thus, the dynamics of these spacetimes can be rigorously described as non-linear, gauge invariant 
and covariant perturbations on a formal FLRW background given by the QL scalars. Since LTB 
spacetimes are compatible with a wide variety of "equations of state" and theoretical assumptions, 
they provide an ideal framework for numerical models of cosmological sources under idealized but 
fully non-linear conditions. As an illustrative example, we briefly examine the formation of a black 
hole in an expanding Chaplygin gas universe. 
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INTRODUCTION. 

The dominant theoretical explanation to account for the cosmic acceleration detected in 
recent observations is based on an elusive source, "dark energy" (DE) that behaves as a 
cosmological constant or as a fluid with negative pressure. While dark matter (DM) 
in galactic halos is assumed to be inhomogeneous (and Newtonian) at the galactic 
scale, dark energy is usually studied by means of FLRW models and/or their linear 
perturbations. The justification is that the dynamical effects of DE are relevant only 
in larger scales (100-300 Mpc), in which the universe appears to be homogeneous or 
nearly so [[T]|. However, as long as the fundamental nature of DM and DE is not known, 
there is no reason to assume a priori that no new valuable information could come from 
studying these sources and their interactions under inhomogeneous, relativistic and fully 
non-linear conditions, at the very least in the galactic and intermediate large scale. 

Since fully general inhomogeneity requires numerical codes of high complexity, we 
offer in this article a compromise by looking at spherically symmetric sources, which are 
obviously more idealized but still useful to examine non-linear phenomena that cannot 
be studied with linear perturbations. The resulting models, "LTB spacetimes", can be 
fully described by autonomous first order evolution equations that can be well handled 
by simple numerical methods. These models are quite general and readily allow for an 
inhomogeneous generalization of a large number of known FLRW solutions. See [|2l[3]| 
for a more comprehensive discussion of the contents of this article. 



LTB SPACETIMES IN THE "FLUID FLOW" DESCRIPTION. 



Spherically symmetric inhomogeneous dust sources are usually described by the well 
known Lemaitre-Tolman-Bondi metric 
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ds^ = -c^dt^ + dr^+R^ (d02 + sin^ 0d(j)^) . (1) 

1 — K 

where R = R{ct,r), R' = dR/dr and K = K{r). A large class of spherically symmetric 
spaectimes follow at once by considering the most general source for ([T]) in a comoving 
frame (w" = 5q), which is the energy-momentum tensor 

r^ = A^MV + p/^''* + ^"^ (2) 

where jU and p are the matter-energy density and the isotropic pressure, h'^^ = u'^u^ + g^^ 
is the induced metric of hypersurfaces ^^^^ orthogonal to u", and H"^ is the symmetric 
traceless tensor of anisotropic pressure. We will call "LTB' spacetimes" to all solutions 
of Einstein's equations for ([T]) and (|2]). 

Besides the scalars ji and p, and the tensor 11"*, the remaining basic covariant objects 
of LTB spacetimes are: 

2,R R' 

= Expansion scalar (3) 

'^M — ^2j^f 7 Ricci scalar of the hypersurfaces (4) 


<yab = ^{aUb)-i^hab Shear tcnsor, (5) 

E"^ = UcUdC"'''''' Electric Weyl tensor (6) 

where R = u"VaR, = /z^V/,, and C"*'''^ is the Weyl tensor. For spherically symmetric 
spacetimes, the symmetric traceless tensors a'^*, IV^^ and can be expressed in terms 
of single scalar functions as 

(j'^b^Y^^ab^ Yl^b^C^'^ab^ ^ab ^ ^ '^ab ^ 

where S"* = h"'^ — 3t]°t]* and r\" = \/lf^5" is the unit vector orthogonal to and to 
the 2-spheres orbits of S0(3) parametrized by (0,^). The field equations G"^ = kT"^ 
(with K = 8kG/c'^) for Q and Q are 

KiiR^R' = [R{R^ + K)]\ (8) 
KpR^R' = -^[R{R^+K) + 2R^R\\ (9) 
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while from ([3]) and (|7]), we obtain for <f and E 
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Bearing in mind (|7]), all covariant objects (scalars and proper tensors) in LTB spacetimes 
can be fully characterized by the local covariant scalars {/i, p, 0, E, ^M}. Given 
the covariant "1+3" slicing afforded by u", their evolution can be completely determined 
by a "fluid flow" description of scalar evolution equations for these scalars (as in 
BmSmi). However, we will consider another covariant scalar representation. 

QUASI-LOCAL (QL) VARLVBLES 

The Misner-Sharp quasi-local mass-energy function, is a well known invariant 
in spherically symmetric spacetimes ifTOl [TTl [T2l [T3l . For LTB spacetimes ([l])-([2]) it 
satisfies the equations 

2^' = K^R^R', (12) 
2^ = -K{p-2^)R^R. (13) 



Comparing (12) with the field equation ([8]) suggest obtaining an integral expression for 



^ that can be related to R and R. This integral along the ^''^^ exists and is bounded 
if the integration domain contains a symmetry center f\3]. Assuming as integration 
domain a spherical comoving region ^ = x C where is the unit 2-sphere, 
^ = {x E R\0 < X < r} and x = marks a symmetry center so that M{ct, 0) = for all 
f, we integrate both sides of ([8]) and also ( [T2] ). This allows us to define a scalar jU* as 

K _ 1J{ _ K foilR^R'dx _R^ + K 

3^*^~R^~3 JoR^R'dx ^ ^ 

where Jq ..dx = J^Zq ..dx. This integral definition of /i*, which is related to /i and to the 
quasi-local mass-energy function, motivates us to define the following: 

Quasi-local (QL) scalar map Let X{^) be the set of all smooth integrable scalar func- 
tions in ^. For every A eX{&), the quasi-local map is defined as 

f'^ /I d2d/j 

^.-.Xm^Xm, ^* = MA)= )rj^2R,^ - (15) 

The scalar functions A* : ^ ^ R that are images of will be denoted by "quasi-local" 
(QL) scalars. In particular, we will call A* the QL dual of A. 



Applying the map ( 15 1 to the scalars and in ^ and Q we obtain 



Applying now ( [15] ) to and p, comparing with ([Sjl-Q, and using ([16]), these two field 
equations transform into 
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0* = 



02 



(/i* + 3j9*). 



(17) 
(18) 



which are identical to the FLRW Friedman and Raychaudhuri equations, but among QL 
scalars. These equations can be further combined to yield identically the FLRW energy 
balance equation: 

/i* = -(/i* + p*) 0*. (19) 

We have found the QL duals for the scalars {/i,p,0, ^^}, with the help of (15) the 
remaining covariant scalars {E, ^ can be expressed as deviations or fluctuations of 
jU, p and with respect to their QL duals: 



while (13) becomes 
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EVOLUTION EQUATIONS FOR THE QUASI-LOCAL SCALARS. 

The scalars A and A* are related by the "relative deviations" or "perturbations" 
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A-A, 



(24) 



which leads to an alternative QL scalar representation {A^, 5^^^ that it is fully equiva- 
lent to the local representation {/i, p, 0, '^^^ E, Hence, LTB spacetimes are fully 
determined by evolution equations for the QL scalars. 

It is straightforward to show from ([15]) that the radial gradients of /i*, p^ and ^ can 
be given in terms of the d functions by 
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(25) 



while ( 18 ) and (19) are evolution equations for /i* and 0^. Hence, the evolution equa- 
tions for 5^^^ and 5^®^ follow from the consistency condition: [Aj] ' = [A*] , apphed 



to (18), (19) and (25) for A* = 0*, /i*. The result is the following set of autonomous 
evolution equations for the QL scalar representation {A*, 5^^^: 
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The constraints associated with these evolution equations are simply the spatial gradients 



(25) , while the Friedman equation (or Hamiltonian constraint) is (17). Notice that the 
constraints ( [25] ) follow directly from differentiating the integral definition (15), so by 
using the QL variables we do not need to solve these constraints in order to integrate 

(26) -(|29]). It is straightforwards to prove (see [3J) that the evolution equations (26)-(29) 
and the constraints (17) and (25) are wholly equivalent to the 1+3 evolution equations 
and constraints for LTB models in the "fluid flow" description of Ellis, Bruni, Dunsbury 
and vanEllst[|7l[8l|9l. 



A NON-LINEAR PERTURBATION SCHEME 



The definition (24) and the evolution equations (26)-(29) suggest that 5^"^^ can be 
rigorously defined as spherical perturbations on a formal FLRW "background" state 
given by the A*. Considering the perturbation formalisms developed by Ellis, Bruni and 
Dunsbury [7l[8l|9l and Bardeen [14J, a perturbation scheme based on the 5^^^ can be 
defined rigorously for the FLRW-LTB case in terms of a suitable map between X and 
X, which are, respectively, the sets of smooth integrable scalar functions in S (FLRW 
model) and S (the "perturbed" lumpy LTB model). For all covariant FLRW scalars A G X 
(we denote FLRW objects with an over-bar) this map is 

^■.X^X, ^{A) = J,{A)=A,eX, 5^^^ = ^^^ (31) 

and characterizes QL scalars (which are LTB objects satisfying FLRW dynamics) as the 
"background model" in LTB spacetimes. Following Dunsbury, Ellis and Bruni EKHl, a 
perturbation scheme on FLRW cosmologies is covariant if S is described by the "1+3" 
fluid flow variables off7^,'8',P|. Although our description of LTB spacetimes is not based 
on these scalars (the local covariant scalars), it is still covariant because /i*, and 0*, 
are themselves covariant scalars by virtue of their connection with the invariants R 
and their derivatives in ( [14] ), ( [T6] ) and ( [23] ) (see UTTll ). Hence, the formalism associated 



with (31 ) is covariant (see [i2ii3]l). 



Also, by virtue of the Stewart- Walker gauge invariance lemma Q, all covariant ob- 
jects in S that would vanish in the background S (a FLRW cosmology in this case) 
are gauge invariant (GI), to all orders, and also in the usual sense (as in [14J). The back- 
ground variables /i*, p*, 0* do not vanish for 5, hence they are "zero order" GI variables 
to all orders. The quantities in LTB spacetimes that vanish for a FLRW cosmology are 
the tensors 11"^, o"^ and E"^, given by ^ in terms of the scalar functions E and ^ 
in (lOHll I. But from (20H22), these functions are basically the fluctuations jU — jU*, 
p — and — 0*. Hence, from (24) and (25), the perturbation variables 5^f^\ 5^p^ and 
5^®\ as well as the gradients /i^, p'^ and 0'^, are all "first order" quantities that are GI 
to all orders. Therefore, LTB spacetimes in the QL scalar representation {A*, 5^^^} are. 
expressible as spherical, non-linear GIC perturbations on a FLRW background. In the 
linear limit these perturbations reduce to spherical perturbations in the long wavelength 
approximation and in the synchronous gauge [[SI. 




FIGURE 1. Hubble scalar and density contrast for a Chaplygin gas LTB model. The figure displays 
(a) the Hubble scalar and (b) the density contrast for a Chaplygin gas overdensity smoothly matched to a 
Chaplygin gas FLRW universe (the matching interface r = /"b = 1 is displayed as a white strip). 




FIGURE 2. Density perturbations. The function S^'-^^ for the Chaplygin gas configuration ofll 



EXAMPLE: THE CHAPLYGIN GAS. 



As in any perturbative approach, we need to impose an "equation of state" (EOS) 
between and /i* to determine the background subsystem ([26])-([27]). Such a choice 



of an EOS also determines the perturbation equations. Evidently, the plausibility of 
the non-linear perturbation formalism must be tested and judged for different EOS 
according to its predictions. We consider as an example the Chaplygin gas as a fully 
relativistic and inhomogeneous LTB spacetime [|2l. This source has been examined 
under inhomogeneous conditions only in a Newtonian context ifTSl . in terms of non- 
linear but relativistic approximations (T6\, or as a static object [17l. 

The LTB spacetime associated with the "standard" Chaplygin gas EOS [[TSl is that 
characterized by the following the following EOS 



P* = -- 



a 



w — — 
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(32) 



where a is a constant. Following (25 1, the relation between pressure and density fluctu- 
ations and perturbations is then 



P-P* 



a 



ill -11* 



(33) 



The evolution equations are simply (26)-([29]) specialized for the EOS (32) and with 
5^P^ given by (ml. Since, from (24[, we have p = p^[l + 5^p^] and /i = ^u^l + 5^^\ 



equations (33) imply the relation 
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[l-(5 



(34) 



which can also be interpreted as a first order virial correction to the FLRW EOS con- 
taining squared fluctuations that convey the effect of long range interactions. This type 
of correction is qualitatively analogous to that arising in the ideal gas under a self- 
gravitating regime in Newtonian systems [|T9l (see also Q). As long as we ignore the 
fundamental physics of the Chaplygin gas, we cannot rule out the possibility that we 
might be describing important non-local effects by using the EOS ([32]) and by having 



the local variables p, ji given by expressions like (34). 



We present here three graphs obtained from the numerical solution of (26)-(29) for 
the Chaplygin gas "top hat" model, constructed by smoothy matching a section of 
a Chaplygin gas LTB model (0 < r < r^) with a Chaplygin gas spatially flat FLRW 
universe (r > r^), with r/, = 1. As shown in panel (a) of figure [T] the QL Hubble scalar 
Mi — 0*/3 passes from infinite values at an initial singularity for all r. In the overdensity 
region we see that ^ — oo, indicating a collapse to a black hole, while in the FLRW 
region ^ tends to a positive constant that can be identified with the A value that FLRW 
Chaplygin models tend to asymptotically. In panel (b) of figure [T] we plot the "density 
contrast" ratio of the QL density /i* to the FLRW density fX. The ratio diverges as inner 
layers of the overdensity collapse into the black hole while external layers blend into 
the cosmic background as r ^ 1. Figure |2] displays the "exact" non-linear perturbation 



d^f^\ This function is negative (overdensity) and close to zero near the center and the 
matching interface (where radial gradients are small). It is significantly different from 
zero in the areas where the radial gradient is large and 5^'^-' ^ — 1 as inner layers collapse 
into a black hole and at the initial singularity for all layers (not shown). A more detailed 
numerical study of the Chaplygin gas LTB model was undertaken in ||2l. 
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